Consequences of Flooding on Spectral Statistics 
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We study spectral statistics in systems with a mixed phase space, in which regions of regular and 
chaotic motion coexist. Increasing their density of states, we observe a transition of the level-spacing 
distribution P(s) from Berry-Robnik to Wigner statistics, although the underlying classical phase- 
space structure and the effective Planck constant h c g remain unchanged. This transition is induced 
by flooding, i.e., the disappearance of regular states due to increasing regular-to-chaotic couplings. 
We account for this effect by a flooding-improved Berry-Robnik distribution, in which an effectively 
reduced size of the regular island enters. To additionally describe power-law level repulsion at small 
spacings, we extend this prediction by explicitly considering the tunneling couplings between regular 
and chaotic states. This results in a flooding- and tunneling-improved Berry-Robnik distribution 
which is in excellent agreement with numerical data. 

PACS numbers: 05.45.Mt, 03.65.Sq 



I. INTRODUCTION 

The universal relation between the statistics of quan- 
tum spectra and classical mechanics is a fundamental 
cornerstone of quantum chaos: For systems with regu- 
lar dynamics it was conjectured that spectral statistics 
show Poissonian behavior In contrast, systems with 
chaotic dynamics should be described by random matrix 
theory [2|, |3|, which can be explained in terms of peri- 
odic orbits [3t9- For generic Hamiltonian systems with 
a mixed phase space, in which disjoint regions of regular 
and chaotic motion coexist, universal spacing statistics 
were obtained by Berry and Robnik j7j]. Their deriva- 
tion is based on the semiclassical eigcnfunction hypothe- 
sis [8j-il0j , which states that eigenfunctions of a quantum 
system scmiclassically localize on those regions in phase 
space a typical that orbit explores in the long-time limit. 
For regular states in one-dimensional systems this corre- 
sponds to the WKB quantization condition [TJ [l2| 



(1) 



p dq = h e g ( m + - 



It shows that the regular state, labeled by the quantum 
number m, localizes on the quantizing torus C m which en- 
closes the area h e s {jn + |) in phase space. On the other 
hand the semiclassical eigenfunction hypothesis implies 
that chaotic states uniformly extend over the chaotic re- 
gion of phase space. Assuming that the disjoint regu- 
lar and chaotic regions give rise to statistically uncorre- 
lated level sequences, one obtains the Berry-Robnik level- 
spacing distribution Q ; see Fig. [T] (dash-dotted lines) . 

The assumption of uncorrelated regular and chaotic 
level sequences does not hold in the presence of dynami- 
cal tunneling [l3l - l26j . which quantum mechanically cou- 
ples regular and chaotic states. If such tunneling cou- 
plings are small, regular eigcnstates will typically have 
tiny chaotic admixtures and vice versa. The influence 



of such weak couplings on spacing statistics can be de- 
scribed perturbatively [27H3l| . Based on this description 
a tunneling-improved Berry-Robnik distribution was de- 




FIG. 1. (Color online) Level-spacing distribution P(s) of the 
model system (see Sec. [TTJ) for h e s » 1/13. The numerical 
data (black histogram) is compared to the flooding-improved 
Berry-Robnik distribution (red dashed lines) , Eq. ([23} , as well 
as to the flooding- and tunneling-improved Berry-Robnik dis- 
tribution (green solid lines), Eq. (|29|l . for system sizes (a) 
M — 1 (weak flooding) and (b) M = 6765 (strong flood- 
ing); M is introduced in Sec. Ill Al For comparison the Wigner 
distribution (dotted lines) and the Berry-Robnik distribution 
(dash-dotted lines) are shown. The insets show averaged 
Husimi functions of chaotic eigenstates. 
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rived recently, which explains the power-law distribution 



of small spacings in mixed systems [31 



For systems with a large density of states, it is observed 
(l5l |32h3t| that a regular WKB state strongly couples to 
many chaotic states. As a consequence, the correspond- 
ing regular eigenstate disappears and chaotic eigenstates 
penetrate into the regular island, ignoring the scmiclassi- 
cal eigenfunction hypothesis. This effect is called flooding 
[33l Hi} . It causes the number iV* of regular eigenstates 
that actually exist in the regular island to be smaller than 
the number A r sc expected from the semiclassical eigen- 
function hypothesis. In Refs. [H, HH it was found that 
in addition to the WKB quantization condition (fTJ) the 
regular state on the mth quantizing torus exists only if 



7m < 



1 



(2) 



Here, 7 m is the tunneling rate, which describes the initial 
exponential decay of the mth WKB state to the chaotic 
region. The Heisenberg time th, c = h e s/A c is the ratio 
of the effective Planck constant h e g and the mean level 
spacing of the chaotic spectrum A c . 

In this paper we study the consequences of flooding on 
spectral statistics in systems with a mixed phase space. 
With increasing density of states we observe a transition 
of the level-spacing distribution from Bcrry-Robnik [see 
Fig- H{a)] to Wigner statistics [see Fig. QJb)], although 
the underlying classical phase-space structure and h e g re- 
main unchanged. This transition is demonstrated quan- 
titatively for model systems with a simple phase-space 
structure, but it is expected to hold for generic sys- 
tems with a mixed phase space. In order to explain 
the transition, we introduce a flooding-improved Bcrry- 
Robnik distribution which takes into account that only 
< iV r sc regular states survive in the regular region. 
We find good agreement with numerical data; see Fig. Q] 
(red dashed lines). We unify this intuitive prediction 
with the tunneling-improved Bcrry-Robnik distribution 
PHI ], which explicitly considers the tunneling couplings 
between regular and chaotic states. This results in a 
tunneling- and flooding-improved Bcrry-Robnik distribu- 
tion, which excellently reproduces the observed transition 
from Bcrry-Robnik to Wigner statistics as well as the 
power-law level repulsion at small spacings; see Fig. Q] 
(green solid lines). 

This paper is organized as follows: In Sec. [TT] we in- 
troduce a family of model systems. Their level-spacing 
distribution is studied in Sec. IIIIl where we demonstrate 
the transition from Berry- Robnik to Wigner statistics nu- 
merically and explain it by the flooding of regular states. 
We conclude with a summary in Sec. IIVI 



II. MODEL SYSTEM 

In this section we introduce a family of model systems 
for which the consequences of flooding can be studied in 
detail. 



A. Classical dynamics 

We consider systems with a mixed phase space where 
classically disjoint regions of regular and chaotic motion 
coexist. As examples we choose one-dimensional kicked 
systems, described by the classical Hamilton function 



H(q,p,t)=T{p) + V(q)J2^-n), 



(3) 



where T(p) is the kinetic energy and the potential V(q) 
is applied once per kicking period. The dynamics of such 
systems is determined by the stroboscopic mapping M. 
of the positions and the momenta (q n ,p n ) at times t = n 
just after each kick [37| . 

M : (q n+1 ,p n+1 ) = (q n + T'(p n ),p n - V'(q n+1 )). (4) 

We design the example systems similar to Refs. [2l|, 
l23l [H, [HI by the piecewise linear functions 



-l + Si(p+l/4) for pe } - 1/2, 0[, 



(5) 



+l-s 2 {p-l/A) for pe]0,l/2[, 
v'(q) = -rq-(l-r)[q + l/2\, (6) 

where \x\ is the floor function and t'(p) is periodically 
extended. Smoothing the functions t'(p) and v'(q) with 
a Gaussian G e (z) = cxp(— z 2 /2e 2 )/v / 27re 2 , one obtains 
analytic functions 



/oo 
dzt'(z)G e (p-z), 
-00 

V'(q) = / dzv'(z)G e (q-z). 



(7) 
(8) 



By construction these functions have the periodicity 
properties 



T'(p + k)=T'(p), 
V'(q + k) = V / (q)-k, 



(9) 
(10) 





FIG. 2. (Color online) Phase-space portrait of the model sys- 
tem, Eq. (j4|. For one unit cell M = 1 (a) the regular island 
(red lines) is embedded in the chaotic sea (blue dots). For 
systems with M > 1 (b) the phase space consists of M such 
unit cells side by side. The arrows indicate the transport in 
the regular islands and in the chaotic sea. 
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for k G Z. This allows to consider the map M on a 
torus, i.e., (q,p) G [-M/2, M/2[x [-1/2, l/2[ with peri- 
odic boundary conditions and M G N. Due to the choice 
of T'(p) and V'(?), the dynamics is equivalent in each 
unit cell of phase space with q E [k — 1/2, fc + l/2[ and 
k E Z; see Fig. [2] In the following we choose the pa- 
rameters si G [5,20], s 2 = 2, r = 0.46, and e = 0.005 
such that each unit cell has a regular island centered 
at (<jk,p) = (fc, 1/4). The area of one such island is 
A r rj 0.32, which equals the relative size of the regular 
region in phase space. 

Since the islands are transporting to the next unit cell 
in the positive q direction, i.e., A4(qk,p) = (qk+i,p), the 
center of each island is a fixed point of the Mth iterate of 
the map, M. M (q~k,P) = (<lk,P)- The surrounding chaotic 
sea has an average drift in the negative q direction as the 
overall transport of the system is zero [111 H(|; see Fig. [2] 
Quantum mechanically this transport suppresses the lo- 
calization of chaotic eigenstates. In our model systems 
the hierarchical regions around the regular islands are 
sufficiently small, and also the effects of partial transport 
barriers and nonlinear resonance chains are irrelevant to 
the numerical studies. 



B. Quantization 

The quantum system is given by the time-evolution 
operator over one period of the driving, 

U = cxp (-J-V{q)\ cxp (--Lt(0) J ; (11) 
V n eS J V "-off / 

see, e.g., Refs. [IH, HI], [42j]. Quantizing the map M 
on a two-torus induces the Bloch phases 9 q and 6 P 
[42I |43| which characterize the quasi-periodicity condi- 
tions on the torus. The Bloch phase 9 p is limited by 
M(6 p + N/2) G Z because of the periodic boundary 
conditions, whereas 9 q G [0, 1[ can be chosen arbitrar- 

ily (Mill- 

Due to the quantization on a compact torus the effec- 
tive Planck constant h c g = 2nh c g is determined by the 
number of unit cells M and the dimension of the Hilbcrt 
space N, 

M 

h cS = ^. (12) 

Here N G N is a free parameter of the quantization and 
the semiclassical limit is reached for h e g — > 0. Note that 
M and N are chosen by continued fractions of h e g = 
l/(d + u) with a = (y/E — l)/2 being the golden mean 
and d G N. This ensures that h e g = M/N is as irrational 
as possible If M and N were commensurate the 

quantum system would effectively be reduced to less than 
M cells. In the following we choose d = 12, leading 
to (M,N) = (1,13), (21, 265), (610, 7697), (6765, 85361), 
such that the effective Planck constant is approximately 
fixed at h c f[ sa 1/13. 



The eigenvalue equation 

J7|0„)=e i ^|0„) (13) 

gives N eigenphases <fr n G [0, 2tt[ with corresponding 
eigenvectors \4> n )- For fixed h e g it is possible to tune 
the density of states by varying AI and N, i.e., for in- 
creasing M, N with approximately constant h e g = M/N 
the density of states rises and flooding becomes more and 
more prominent, as will be discussed in Sec. 1111 Bl In or- 
der to numerically solve the eigenvalue equation (|13p for 
N > 10 4 we use a band-matrix algorithm, see the Ap- 
pendix. 

III. SPECTRAL STATISTICS AND FLOODING 

In this section we study the consequences of flood- 
ing on spectral statistics. In Sec. IIII AI we consider the 
model systems introduced in Sec. [ill Increasing their den- 
sity of states (M —> 00) at fixed h e g gives the flooding 
limit for which we obtain a transition of the level-spacing 
distribution P(s) from Berry-Robnik to Wigner statis- 
tics. In Sec. IIII Bl we discuss flooding of regular states. 
Based on this discussion, we introduce the flooding- 
improved Berry-Robnik distribution Pg(s) in Sec. IIII Cl 
which intuitively explains how the flooding of regular 
states causes the transition from Berry-Robnik to Wigner 
statistics. In Sec. IIII Dl we unify this prediction with 
the results of Ref. |3l|. leading to the more sophisticated 
flooding- and tunncling-improved Berry-Robnik distribu- 
tion Pfti(s). This distribution additionally accounts for 
the effects of level repulsion between regular and chaotic 
states. In Sec. IIII El we consider three limiting cases in 
which level repulsion vanishes. In particular we discuss 
that the semiclassical limit, h c g —> with fixed M, leads 
to the standard Berry-Robnik statistics, while Wigner 
statistics are obtained in the flooding limit considered in 
this paper. 

A. Spacing statistics of the model system 

We investigate the spectral statistics of the model 
systems introduced in Sec. [II] numerically. In order 
to increase the statistical significance of the spectral 
data, we perform ensemble averages by varying the pa- 
rameter si of the map; see Eq. ([5J). This modifies 
the chaotic dynamics but leaves the dynamics of the 
regular region unchanged. Also the Bloch phase 6 q 
is used for ensemble averaging. For the parameters 
(M,N) = (1,13), (21, 265), (610, 7697), (6765, 85361) we 
choose 50, 50, 10, 4 equidistant values of si in [5, 20] and 
400, 19, 10, 1 equidistant values of 6 q in [0, 1[, respectively. 
For each choice the ordered eigenphases </> n give the un- 
folded level spacings 

TV 

S n ■= —{4>n+l - 4>n)- (14) 
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Assuming an uncorrclated superposition of regular and 
chaotic subspectra corresponding to disjoint regular and 
chaotic regions in phase space, these spacings are ex- 
pected to follow the Berry- Robnik distribution The 
only relevant parameter of this distribution is the density 
of regular states which semiclassically equals the relative 
size of the regular region in phase space, A r ; see Eq. (TT51) . 
This gives the standard Bcrry-Robnik distribution 



Pbr(s) 



ds 2 



cxp(-^ r s) erfc ^—(1 - A r )s 



(15) 



For purely chaotic systems one has A r = such that the 
Wigncr distribution P c (s) = (7rs/2)e _7rs l A is recovered. 
For purely regular systems one has A r = 1, giving the 
Poisson distribution P r (s) = e _s . 

For the model systems introduced in Sec. [TT| one has 
A T S3 0.32, such that Eq. (fTS")) predicts the same level- 
spacing distribution for all system sizes M. This is in 
contrast to our numerical findings, which show a transi- 
tion of the level-spacing distribution from Bcrry-Robnik 
to Wigner statistics with increasing system size M and 
fixed h c ff. In Fig. [3] numerical results for the level-spacing 
distribution of the model systems are shown as black his- 
tograms. For the case of only one unit cell [Fig. EJa)] 
the level-spacing distribution roughly follows the Bcrry- 
Robnik distribution (dash-dotted line). With increase 
of the system size to M = 21 unit cells [Fig. Elb)] the 
level-spacing distribution shows global deviations from 
the Bcrry-Robnik distribution. For even larger system 
sizes [Figs. [S^c) and EJd)] wc observe a transition to 
the Wigner distribution (dotted line). This transition 
is caused by flooding of regular states, which we discuss 
in the following section. 



B. Flooding of regular states 

We now show how the number of regular and chaotic 
states is modified in the presence of flooding. According 
to Eq. (TTJ) each regular state occupies an area /i c ff hi 
phase space. Hence, the maximal number of quantizing 
tori m max per island is given by 



■;?).„ 



(16) 



and the quantum number m runs from to TO max — 1- 
Since we consider a chain of M islands there are M reg- 
ular levels supported by the mth quantizing tori of the 
M islands. Hence, we semiclassically expect 



A; s 



C M M A r M/h e s 



(17) 



regular states supported by the M regular islands of size 
A r . The semiclassically expected density of regular states 
p s r c is therefore given by the relative size of the regular 
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FIG. 3. (Color online) Level-spacing distribution P(s) of the 
model system for h e g ~ 1/13. The numerical data (black 
histograms) show a transition from the Berry-Robnik distri- 
bution (dash-dotted lines) to the Wigner distribution (dotted 
lines) with increasing system size M = 1, 21, 610, and 6765 
(a)-(d). These data are compared to the flooding-improved 
Berry-Robnik distribution (red dashed lines), Eq. 1)230 . as well 
as to the flooding- and tunncling-improved Berry-Robnik dis- 
tribution (green solid lines), Eq. (|29l) . The insets show the 
same distributions on a double logarithmic scale. 



region, 



N sc 

S3 A*. 

N 



(18) 

chaotic states and 



Similarly we expect A c sc = N - N? 
P? = 1 - P?- 

Due to dynamical tunneling, regular and chaotic states 
are coupled. The average coupling of the regular states 
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localizing on the mth quantizing tori to the chaotic states states, Eq. (fl"8|) . and using Eqs. (fT2|) and (ITTl) we obtain 



is given by the typical coupling v m [31| . It is determined 
by the tunneling rate j m which describes the initial ex- 
ponential decay of the mth regular WKB state to the 
chaotic sea, 



N 
2^ 



A 7 ; 



l 

2^ 



In 



Ksp s c 



M. 



(19) 



We compute the system specific tunneling rates j m nu- 
merically [23|. They depend only on Planck's constant 
h e ff and the classical phase-space structure of one regu- 
lar island, which are fixed in our investigations. Hence, 
the factor \/j m / (h spl c ) in Eq. (Q~9|) is constant for our 
model systems and the typical coupling v m is tunable by 
the system size M. Note that the couplings v used in 
Refs. [HI, HH differ by the factor p^ c from our definition, 
Eq. (p~9]> . due to a different choice of dimensionless units. 

In Ref. [HI it was shown that in addition to the WKB 
quantization condition ([T]) regular states exist on the 
?nth quantizing tori only if the tunneling rate 7 m is 
smaller than the inverse Heisenberg time of the chaotic 
subsystem, 7 m < 1/th, c , Eq. @. Using Eq. ([TO]) and 
th,c = h c ff/A c = N^ c we rewrite this existence criterion 
in terms of the typical coupling, 



< 



1 



2tt^ 



(20) 



If the existence criterion (|20|) is fulfilled, the typical cou- 
pling of the WKB states on the mth quantizing tori is 
smaller than the chaotic mean level spacing and the cor- 
responding regular cigenstates exist. If v m increases be- 
yond this threshold, the regular states on the mth quan- 
tizing tori effectively couple to an increasing number of 
spectrally close chaotic states. Consequently the corre- 
sponding regular eigenstates disappear. This process is 
called flooding of regular states [H, H3, 0, EH . Thus 
for large typical couplings v m the number iV r of regu- 
lar states which actually exist in the regular islands is 
smaller than the semiclassically expected number N* c of 
regular states. The quantizing tori of the iV r sc — N* reg- 
ular states which violate Eq. (|2U)) are flooded by chaotic 
states in phase space. Note that for our model systems 
the relation vq < V\ < v-2 < . . . holds, such that the quan- 
tizing tori are flooded in the order of decreasing quantum 
number m from the border to the center of the regular 
islands. 



Pi 



(21) 



771=0 



This expression shows that each quantizing torus semi- 
classically contributes one Planck cell to the density of 
states. To compute the density of regular states p[ in the 
presence of flooding we include only those quantizing tori 
in the sum in Eq. pip for which the existence criterion 
© holds, and obtain [H, |H| 



f 

Pi 



771=0 



'ciT 



i - e [v. 



2rrpg 



(22) 



In Fig. |4] the density of regular states p\ is shown for 
the example system at h e tj ~ 1/13 versus the system size 
M on a logarithmic abscissa. It decreases with increas- 
ing system size M and has a step whenever a typical 
coupling v m equals l/(2npl c ). The semiclassical density 
of regular states p^ c is an upper limit. In the spirit of 
Eq. (fTS]) , w A r , we interpret the density of regular 
states p\ in the presence of flooding by an area in phase 
space. For M = 1, 21, 610, and 6765 the insets of Fig. H 
show this area enclosed by a gray torus. In addition the 
averaged Husimi functions (bright /yellow to darker/red 
color scale) illustrate that the surviving regular states are 
localized in this area p\, which decreases in the flooding 
limit M — > oo. Already for the system with one unit cell, 
71/ = 1, we find that p{ is smaller than its semiclassical 
expectation because the outermost regular state of 
quantum number m = 3 violates the existence criterion 




C. Flooding-improved Berry- Robnik distribution 

We now introduce a flooding-improved Bcrry-Robnik 
distribution which takes the flooding of regular states 
into account. For that purpose we compute the density 
of regular states p\ in the presence of flooding. Start- 
ing from the semiclassically expected density of regular 



FIG. 4. (Color online) Density of regular states p\ in the 
presence of flooding, Eq. ()22|) . vs. the system size M for 
h e fi ~ 1/13 on a logarithmic abscissa. For comparison, the 
semiclassically expected density of regular states p° c is shown 
(green dotted line). The insets illustrate the classical phase 
space where the gray tori enclose the area pi for M — 1, 21, 
610, and 6765. In addition, the averaged Husimi function of 
all chaotic eigenstates folded into the first unit cell is shown. 
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Note that the amount by which a regular state is 
flooded can also be described by smooth functions, e.g., 
by the fraction of regular states [HI or the asymptotic 
flooding weight [45j . However, they do not provide a 
significant advantage for our investigations. 

To obtain a description of the level-spacing distribution 
which includes flooding, we now use prediction (|22|) for 
the density of regular states p[ instead of the relative 
size of the regular region A r as the relevant parameter in 
Eq. ([15]). With p{ := 1 - pi this leads to the flooding- 
improved Berry-Robnik distribution 



-Pfi(s) = ^ jexp(-^s) erfc (^j-pls 



(23) 



which is our first main result. In Fig. [3] we com- 
pare the numerically determined nearest-neighbor level- 
spacing distribution to the analytical prediction (|23[) 
for our model system. With increasing system size M 
and fixed h e g we find a global transition of the level- 
spacing distribution from Berry-Robnik to Wigner statis- 
tics. This global transition is well described by the 
flooding-improved Berry-Robnik distribution, Eq. (|23|) . 
It is a consequence of flooding, which reduces the den- 
sity of regular states below its semiclassical expectation, 
pl < j°r C - According to Eq. (|TS)) this can be interpreted 
as a flooding induced decrease of the regular region in 
phase space. In the limit M — » oo the regular islands are 
completely flooded and no regular state exists. Hence, 
the Wigner distribution is obtained. 

Note that even for the case of only one unit cell [see 
Fig. EUa)] non-zero couplings v m exist such that the 
numerical data are better described by the flooding- 
improved Berry-Robnik distribution, Eq. (|23[) . than by 
Eq. CED. 

At small spacings deviations between numerical data 
and the flooding-improved Berry-Robnik distribution are 
visible. They occur due to level repulsion between the 
surviving regular and the chaotic levels, which is not con- 
sidered within this approach and will be incorporated in 
the following section. 



D. Flooding- and tunneling-improved 
Berry-Robnik distribution 

We now unify the flooding-improved Berry-Robnik dis- 
tribution, Eq. (|2"3"|) . with the tunneling-improved Berry- 
Robnik distribution [3l|. The resulting flooding- and 
tunneling-improved Berry-Robnik distribution allows us 
to describe both the effect of flooding and the power-law 
level repulsion at small spacings. The derivation is done 
along the lines of Ref. 31] . We incorporate the effects 
of flooding into this theory by replacing the number of 
regular states 7V r sc with the number of surviving regular 
states Nf. which fulfill the existence criterion (|20l) . The 
other regular states, which fulfill the WKB quantization 
condition ([T]) but have strong couplings to the chaotic 



sea, v m > l/(27rpp C ), are assigned to the chaotic subspec- 
trum. Level repulsion is then modeled by accounting for 
the small tunneling couplings v m between the N* surviv- 
ing regular states and the chaotic states perturbatively. 

Following Refs. @, [H, HI EH the flooding- and 
tunneling-improved Berry-Robnik distribution Pm{s) 
consists of three distinct parts: 



-Pfti(s) =p r -r(s) +Pc-c(s) + p r - c {s). 



(24) 



Here p r _ r (s) describes the contribution of level spacings 
between two regular levels, p c _ c (s) the contribution of 
level spacings between two chaotic levels, and p r - c (s) 
the contribution of level spacings formed by one regu- 
lar and one chaotic level in the superposed spectrum. In 
our model systems the number of quantizing tori TO max is 
small, e.g., m max w 4, and the M regular levels with the 
same quantum number m are eq uisp aced with distance 
N/M in the unfolded spectrum [46[. Hence, the regu- 
lar levels do not follow the generic Poissonian behavior 
occurring for large m max , but are well separated, 



Furthermore, 



p r _ r (s) w 0. 



p c - c (s)=P c {s)[l-p?sl 



(25) 



(26) 



where P c (s) is the Wigner distribution, which describes 
the probability of finding a spacing s in the chaotic sub- 
spectrum. The second factor [1— p®°s] describes the prob- 
ability of finding a gap in the regular subspectrum. For 
the third term in Eq. (|24|) one finds (3~0 | 



Nl-l 



Pr-^)=P^ c (s)^t^ X 



1 771— 



S 
2Vn 



(27) 



with X(x) := y/n/2x exp(— x 2 /4)/ (a: 2 /4), where Iq is 
the zeroth-ordcr modified Besscl function and v m = 
v m /\/l — 2ir(p* c v m ) 2 . The contribution of the zeroth- 

order regular-chaotic spacings, Pr-c( s )> ^ s given by 



p^ c ( S ) = 2p^ c exp 



API 



(28) 



Using Eqs. (|2"5]l. (gg]) , and (gTJ) in Eq. flU}, we obtain the 
flooding and tunneling improved Berry-Robnik distribu- 
tion 



P fti ( s) = P c(s )[l_ /9fs]+J ,(o_) c(s) 1 j2 v ^x(J- 

1 m— v 

(29) 

which is our final result. In Eq. (|29p one has to sum over 
the iV r < A^ c regular states, which fulfill the existence 
criterion (f20)) . This selection of the regular states takes 
flooding into account. In addition power-law level repul- 
sion at small spacings is described by the last term of 
Eq. (EU). 
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In Fig. [3] we compare the numerical results for the 
level-spacing distribution to the flooding- and tunneling- 
improvcd Bcrry-Robnik distribution, Eq. (|29l) (green 
solid lines). We find excellent agreement. The global 
transition of the level-spacing distribution from the 
Berry- Robnik distribution in Fig. EJa) for a system with 
one unit cell to the Wigner distribution in Fig. |3Id) for 
a system with M = 6765 is well described. This transi- 
tion is caused by the disappearance of regular states due 
to flooding. Furthermore, the flooding- and tunneling- 
improved Berry- Robnik distribution, Eq. (|29[) . repro- 
duces the power-law level repulsion of P(s) at small spac- 
ings, which is caused by small tunneling splittings be- 
tween the surviving regular and chaotic states. This can 
be seen particularly well in the double logarithmic insets 

of Fig. m 



E. Limiting cases 

Depending on the interplay between the flooding limit 
M — > oo and the semiclassical limit h e s — > 0, we identify 
three cases in which the tunneling corrections of Sec. MI Dl 
are insignificant. 

Case (i) is the flooding limit with fixed h e g and M —> 
oo in which all regular states are flooded. Asymptotically 
one obtains the Wigner distribution; see Fig.(3Jd). Note 
that a further increase of the system size after all regu- 
lar states have been flooded completely, may lead to the 
localization of chaotic eigenstates which affects spectral 
statistics [47l.l48j. 

Case (ii) considers the semiclassical limit h e s — > for 
fixed system sizes M. In this case both flooding and tun- 
neling corrections vanish due to exponentially decreasing 
tunneling couplings. Hence, the spacing statistics tend 
toward the standard Berry-Robnik distribution, Eq. (fTS")) 

and the 



oo arc coupled such that p\ is con- 



In case (iii) the semiclassical limit h e g 
flooding limit M 
stant and smaller than p^ c . In this limit flooding is 
present yet the tunneling corrections at small spacings 
vanish. In this case the spacing statistics are given by the 
flooding- improved Berry-Robnik distribution, Eq. (j2"5)l . 

In Fig. [5] we illustrate spectral statistics in the limit of 
case (iii) . We consider the model systems for (M, N) — 
(5,63), (34,735), (89,2458), (1597,77643) such that the 
density of regular states is fixed, p\ « 0.58p^ c , and 
fteff « 1/13,1/22,1/28,1/49 decreases. Both the nu- 
merical data and the flooding- and tunneling-improved 
Berry-Robnik distribution, Eq. (f2"9")l . tend toward the 
flooding-improved Berry-Robnik distribution, Eq. (f2l?|) . 
The vanishing influence of the tunneling corrections at 
small spacings is particularly visible in the insets, which 
show the spacing distributions on a double logarithmic 
scale. 
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FIG. 5. (Color online) Level-spacing distribution P(s) of 
the model system at fixed density of regular states pi ~ 
0.58/5?°. The numerical data (black histograms) are com- 
pared to the flooding-improved Berry-Robnik distribution 
(red dashed lines), Eq. (|23p . as well as to the flooding- and 
tunneling-improved Berry-Robnik distribution (green solid 
lines), Eq. {29]), for (M, N) = (5,63), (34,735), (89,2458), 
and (1597, 77643) (a)-(d), corresponding to h cB « 1/13, 1/22, 
1/28, and 1/49, respectively. For comparison, the Wigner 
distribution (dotted lines) and the Berry-Robnik distribution 
(dash-dotted lines) are shown. The insets show the same dis- 
tributions on a double logarithmic scale. 



IV. SUMMARY 

In this paper we study the impact of flooding on the 
level-spacing distribution P(s) for systems with a mixed 
phase space. Numerically we find a transition from 
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Berry-Robnik to Wigncr statistics with increasing den- 
sity of states and fixed h e g. We explain this transition 
by the flooding of regular islands. It reduces the density 
of regular states p\ below its semiclassical expectation 
p a Y c which can be interpreted as a flooding induced de- 
crease of the regular region in phase space. Taking this 
into account we derive a flooding-improvcd Berry-Robnik 
distribution, which reproduces the observed transition 
of the level-spacing statistics. We unify this prediction 
with the tunneling-improved Berry-Robnik distribution 
[3l| which includes power-law level repulsion. This gives 
the flooding- and tunneling-improved Berry-Robnik dis- 
tribution, which shows excellent agreement with numer- 
ical data. 

In order to demonstrate the effect of flooding on spec- 
tral statistics, we investigated model systems with a 
simple phase-space structure. However, we expect that 
flooding has similar consequences for systems with a 
generic phase space, which may contain several regu- 
lar and chaotic regions as well as nonlinear resonances, 
larger hierarchical regions, and partial transport barri- 
ers. This expectation is based on the fact that even for 
generic systems tunneling couplings between classically 
disjoint regions exist. Hence, increasing the density of 
states for fixed h e g should still lead to a transition of the 
level-spacing distribution from Berry-Robnik to Wigncr 
statistics. 
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APPENDIX: PERIODIC BAND MATRIX 

Our aim is to calculate the eigenphases <p n from 
Eq. (fT3|) numerically up to Hilbert-space dimension TV ss 
10 5 . This is possible due to a band- matrix algorithm (Hoj 
which was used in Rcfs. [HI, [34| and is presented in the 
following. 

We start with the matrix representation of the sym- 
metrized time-evolution operator, 

j/sym = e -^T(P) c -n^V(q) e "^T(^ ^ 

in the basis of the discretized position states \qk), 

uz m --={q k \u sym h), (3i) 

with q k = h eS (k + p - §) and k, I = 0, 1, . . . , N - 1. This 
matrix has dominant contributions around the diagonal 
and in the upper right and lower left corners, i.e., U sym 
can be approximated by a periodic band matrix. For 



our model systems the width of the band depends on the 
extrema of V'(q). The essential step for computing the 
eigenvalues of U sym is to find a similarity transformation 
from this periodic band matrix to a band matrix. For 
the Hermitian case a similar idea was used in Ref . [5l| . 

Since the kicking potential is symmetric about q = 
for M(9 p + TV/2) G Z, V(qi) = V(q N -i), wc find 

U k,l — U N-l,N-k- 

Hence, the set of eigenvectors \<fi n ), for which we choose 
the phase such that (qo\4> n ) = (<Zo|0n)*, satisfies 

(qi\<t>n) = (qN-l\<f>n)* , (33) 

where the star denotes the complex conjugation and I 
runs from 1 to N — 1. Based on these relations it is 
possible to find a unitary transformation A to a set of 



purely real vectors \ip n ) '■= A | (/>„), given by 

MVO - (qo\4> n ), (34) 

(<72fc-l|-0n) = -j={{q k \<t>n) + (q N -k\(/>n)), (35) 

(<72fc|-0n) = — 7={{q k \<p n ) - {q N -k\(f>n)), (36) 

fev-ilVVi) = (?jv/2 I </!>«)• (37) 



Here, k runs from 1 to (N — l)/2 for odd N or from 
I to (N — 2)/2 for even N and the last row has to be 
considered only for even N. 

We now define a new operator W, given by the unitary 
transformation of U s y nl with A, 

W :=AU aym A- x . (38) 

Its matrix representation W in the basis of position states 
\q k ) has a banded structure with twice the bandwidth of 
the matrix U sym but without components in the upper 
right and lower left corners. Furthermore it is symmetric, 

W k ,i = W l<k , (39) 

with complex matrix elements Wk,i- The unitary trans- 
formation (|38p leads to a new eigenvalue problem with 
the same eigenphases as in Eq. ([TBI . 

W\*p n ) = e i ^\ij n ). (40) 

Numerical standard libraries provide methods for the 
eigenvalue computation of only real symmetric or com- 
plex Hermitian band matrices but not for unitary band 
matrices such as W. Hence, we first calculate the real 
part of the eigenvalues, following from Rc{VK} \ip n ) = 
cos 4> n \tp n } , and afterwards the imaginary part from 
Im{VF} \ip m ) = smcj) m \ip m ). This is possible since the 
eigenvectors \tp n ) are purely real. From these results 
one can recover the eigenphases <f> n by the requirement 
cos 2 <j> n + sin 2 4> m = I . The corresponding eigenfunctions 
\ip n ) can be obtained from Eq. P0"|) by the method of in- 
verse iteration using an LU decomposition. By the map- 
ping of the original eigenvalue problem Eq. (| X3[) to the 
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band matrix form Eq. (|40[) , it is possible to compute both 
eigenvalues and eigenfunctions with a numerical effort 



scaling as N 2 in contrast to the standard diagonalization 
procedures which scale as iV 3 . 
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